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Abstract 

We calculate the frequency of collective modes of a one-dimensional repulsively interacting Fermi 
gas with high-spin symmetry confined in harmonic traps at zero temperature. This is a system 
realizable with fermionic alkaline-earth-metal atoms such as which displays an exact SU(k) 

spin symmetry with n ^ 2 and behaves like a spinless interacting Bose gas in the limit of infinite 
spin components k ^ oo, namely high-spin bosonization. We solve the homogeneous equation 
of state of the high-spin Fermi system by using Bethe ansatz technique and obtain the density 
distribution in harmonic traps based on local density approximation. The frequency of collective 
modes is calculated by exactly solving the zero-temperature hydrodynamic equation. In the limit 
of large number of spin-components, we show that the mode frequency of the system approaches to 
that of a one-dimensional spinless interacting Bose gas, as a result of high-spin bosonization. Our 
prediction of collective modes is in excellent agreement with a very recent measurement for a Fermi 
gas of ^^^Yb atoms with tunable spin confined in a two-dimensional tight optical lattice. 
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I. INTRODUCTION 


Ultracold atomic gases appear to be a versatile tool for discovering new phenomena and 
exploring new horizons in diverse branches of physics. To large extent, this is dne to their 
unprecedented controllability and purity. A vast range of interactions, geometries and di¬ 


mensions is possible: using the tool of Feshbach resonances {l| and applying a magnetic held 
at the right strength, one can control very accurately the interactions between atoms, from 
arbitrarily weak to arbitrarily strong. By using the technique of optical lattices that trap 
atoms in crystal-like structures Q, one can create artihcial one- or two-dimensional envi¬ 
ronments to explore how physics changes with dimensionality. Most recently, it is also able 
to control the number of spin-components: degenerate Fermi gas with high-spin symmetry 
has been observed in alkaline-earth-metal atoms Yb 3|, 1^. 

The ytterbium (Yb) atom has a unique advantage in studying high-spin physics. It has 
a closed-shell electronic structure in the ground state, and hence its total spin is determined 
entirely by the nuclear spin, I. For the fermionic species ^^^Yb, the nuclear spin is / = b/2 
and the atom can be in six different internal states k = 21 + 1 = 6. This gives rise 
to a unique feature of ^^^Yb atom, that is, one simple internal-state-independent s-wave 
scattering length due to the absence of electronic spin in the atomic ground state joj. Thus, 
the system exhibits SU(6) symmetry 7H9j. In such a high-spin systenn possible novel ground 
states and topological excitations have been addressed theoretically THlOj. One dimensional 
(ID) repulsively interacting fermions with sufficiently high spin also behave like spinless 
interacting Bose atoms {ill, a phenomenon that we may refer to as high-spin bosonization. 
Physically, this phenomenon may also occur in two or three dimensions. 

Experimentally, Fermi degeneracy of ^^^Yb and ^^^Yb atoms has been demonstrated 
In particular, in a very recent experiment performed at European Laboratory for 
Non-Linear Spectroscopy (LENS), a Fermi gas of ^^^Yb atoms has been created in one¬ 
dimensional harmonic traps by using optical lattices and its momentum distribution and 


B. 


breathing mode oscillation have been measured 

n 

In this work, motivated by the recent measurement at LENS [3, we investigate a ID 
high-spin Fermi gas with strongly repulsive interactions satisfying SUfft:) symmetry, by using 
exact Bethe ansatz technique beyond the mean-held framework [12|. First, we solve the exact 
ground state of a homogeneous Fermi gas at zero temperature based on Bethe ansatz. To 
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FIG. 1: (Color online) Breathing mode frequency as a function of the number of components at 
the dimensionless interaction parameter Na\j^/a\^ = 0.44 ± 0.08. The solid circles with error bars 
are the experimental data reported by the LENS team [^. The empty squares are the theoretical 
results, with error bars counting for the experimental uncertainty for interaction parameter. The 
thin horizontal line at the right part of the figure shows the theoretical prediction at the infinitely 


large number of components. Source: Adapted from Ref. 
Group. 


0 |. 
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make contact with the experiment, we then consider an inhomogeneous Fermi cloud under 
harmonic conhnement, within the framework of local density approximation (LDA). The 
equation of state of the system and the density distribution are calculated. By solving the 
zero-temperature hydrodynamic equation, we predict the frequency of low-lying collective 
modes. The main outcome of our research is summarized in Fig. [H which shows the predicted 
breathing mode frequency in comparison with the recent data reported by LENS group, for 
a trapped ^^^Yb gas at a dimensionless interaction parameter Na\jj/a\^ = 0.44 ± 0.08 j^. 
Here, N is the total number of atoms in the Fermi cloud, am and aho are the ID s-wave 
scattering length and the length of harmonic oscillator, respectively. We hnd an excellent 
agreement between our theoretical prediction and the experimental data, with a relative 
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discrepancy at a few percents. Fig. 1 has been pnblished in Ref. The pnrpose of this 
paper is to present the details of onr calcnlations, focnsing particnlarly on the high-spin 
bosonization. 

We note that in earlier works ID ninlti-component Fermi gases with large-spin and attrac¬ 
tive interactions have been discnssed in detail {^, 141. However, for repnlsive interactions, 
only homogeneons Fermi gas in weak and strong conpling limits has been considered 15| . 


Comparing with the case with attractive interactions there are no mnlti-component bonnd 
clnsters in repnlsively interacting Fermi gases. 

The paper is organized as follows. In the following section, we describe briefly the model 
Hamiltonian. In Sec. HI, we present the exact Bethe ansatz solntion and discnss the eqnation 
of state and sonnd velocity of a nniform Fermi gas at zero temperatnre. Then, by nsing LDA 
in Sec. IV we determine the density distribntion in the trapped environment. In Sec. V, 
we describe the dynamics of trapped Fermi gases in terms of ID hydrodynamic eqnation. 
The behavior of low-lying collective modes is obtained and discnssed. Finally, we conclnde 
in Sec. VI. 


II. MODEL HAMILTONIAN 


We consider a ID mnlti-component Fermi gas with psendo-spin S' = (k — 1) /2, where k 
(> 2) is the nnmber of components. The fermions in different spin states repnlsively interact 
with each other via the same short-range potential giD6{x). The hrst-qnantized Hamiltonian 
with a total nnmber of atoms N = Yl'i=i (where Ni is the nnmber of fermions in the 
psendo-spin state 1) for the system is 



( 1 ) 


where V = mu‘^x‘^/2 is the harmonic trapping potential for the atom i and u is the trapping 
freqnency. In general, snch a ID Fermi system is created by loading a 3D clond into a tight 
2D optical lattice and separating it into a nnmber of highly elongated tnbes. In each tnbe, 
it is convenient to express ID conpling constant giu in terms of an effective ID scattering 
length. 


giD 


2h^ 

maiD ’ 


( 2 ) 
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FIG. 2: (Color online) Dependence of the quasi-momentum distribution g (x) on the interaction 
strength (a) and on the number of components (b). At x = 0, g{x) increases with decreasing the 
interaction strength or increasing the number of spin-components. While at large x, g{x) saturates 
to l/27r. 


where the effec 
the relation 16 


ive ID scattering length air, is related to the 3D scattering length a^r, by 

m , 


aiD = {I-A— 1 < 0 . 

CtSD 


( 3 ) 


Here, Op = \/h/(mup) is the characteristic oscillator length with transverse frequency Up 
determined by optical lattice depth and 


A=- ^ ’’ ~ 1.0326 
\/2 


( 4 ) 


is a constant. 


III. HOMOGENEOUS EQUATION OF STATE 

Let us hrst address a ID uniform multi-component Fermi gas with symmetric inter¬ 
component interactions, i.e., there is no trapping potential term Vi in the Hamiltonian 
(fTj). In this case, the model Hamiltonian is exactly soluble via the Bethe ansatz technique 
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12|. Focusing on the experiment at LENS, we assume that each component has the same 
number of particles, i.e., Ni = N/k for I = 1,2, ■■■, k = 2S + 1. In free space, we measure 
the interactions by a dimensionless coupling constant 

- ^9id 2 

7 = —^ =-, (5) 

n^^n. —I naiD 

where n is the linear total number density The ground state energy Ehom in the 




thermodynamic limit is given by 

F fc 2 


m 


nyy 


x^g (x) dx. 


( 6 ) 


where 


A = 27 / g (x) dx 


( 7 ) 


and the quasi-momentum distribution g (x) with x > 0 is determined by an integral equation 


9i^) = / 9 {x')dx' 


1=1 

IX 


IX 


.(/A)^ -I- (x — x')^ 


(/A)^ (x -I- x')^ 


( 8 ) 

The above Bethe ansatz equations are very similar to those for attractive interactions Q,Q 
19l. |20| . but without the contribution from K-body bound cluster states. In Fig. [21 we show 
the quasi-momentum distribution g (x) for typical interaction parameters and the number 
of spin components, obtained by solving the integral equation Eq. ([ 8 ]), together with Eq. 

( 0 . 


Once we obtain the ground state energy by using the Bethe ansatz technique, we calculate 
the chemical potential by /ihom = dEy^om/dN and the corresponding sound velocity by Chom = 
\/n{dfihom/dn)/m. For numerical purposes, it is convenient to rewrite these thermodynamic 
quantities in a dimensionless form that depends on the coupling constant 7 only. 


T^hom 
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/^hom 


Chom 
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e(7), 
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2m 
hn 

cicf). 
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( 9 ) 
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( 10 ) 
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FIG. 3: (Color online) Dependence of the uniform gronnd state energy per particle (a), the chemical 
potential (b), and the velocity of sonnd (c) on the dimensionless conpling constant 7 , at several 
number of species as indicated. The energy per particle and sound velocity are in units of the 
energy 7 r^[^^n^/( 2 m)] and the velocity 7 r(^n/m), respectively. Thin solid lines in (a) and (b) are 
the results of a ID spinless interacting Bose gas with the same density n and coupling constant 7 . 
Note that in (b) for chemical potential, 7 is shown in the logarithmic scale. Thus, the saturation 
to the strongly interacting Tonk-Girardeau limit is not obvious. 


These dimensionless functions are related by. 


Qp 

/r(7) = 3e(7)-7^, 

(12) 

/X / / X 7^/^ ( 7 ) 

c(7) = hM7) 2 a, ■ 

(13) 
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It is easy to see that for an ideal, non-interacting mnlti-component Fermi gas, 


T^ideal 

■^hom 


2m 


(ttVSk^) , 




and 


ideal _ 

^hom 


hn 


m 


t/k). 


( 14 ) 

( 16 ) 

( 16 ) 


By nnmerically solving the integral Eqs. (I6])-(IH!), we obtain the gronnd state energy per 
particle, chemical potential, and sonnd velocity as a fnnction of the dimensionless conpling 
constant 7 , as shown in Fig. [3l With increasing the conpling constant, the energy, chemical 
potential and sonnd velocity increase rapidly from the ideal gas resnlts and hnally satnrate 
to the strongly interacting Tonk-Girardean gas limit, as one may anticipate. With increasing 
the nnmber of components k, these thermodynamic qnantities decrease instead. It is inter¬ 
esting that for a snfhciently large nnmber of components, they approach to the resnlts of a 
ID repnlsively interacting spinless Bose gas with the same total density n and the conpling 
constant 7 = —mgiD/{f^n), which are obtained by solving the following integral eqnation 


lal, 




+1 


A2 


B 


-1 


+ (fc - k') 


f{k')dk\ 


( 17 ) 


together with the normalization condition 


Ab ^ 27 / / [k) dk 


( 18 ) 


and the expression for the energy 


t^B 

■^hom 


m 


77-7 

Ab 


k^f (k) dk. 


( 19 ) 


The eqnivalence between high-spin Fermi gas and spinless Bose gas, namely high-spin 
bosonization, has been analytically shown by Yang and Yon ll|. This is an interesting 
connterpart of the ID effective fermionization for strongly interacting particles in one di¬ 
mension 21 1 . 
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FIG. 4: (Color online) Density distributions of a ID trapped multi-component Fermi cloud at 
three interaction parameters = 0.1 (black solid line), 1 (red dashed line), and 10 (blue 

dotted line). The linear density and the coordinate are in units of the peak density and 

Thomas-Fermi radius xtf,k of an ideal Fermi gas, respectively. 
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IV. DENSITY DISTRIBUTION IN HARMONIC TRAPS 


li, 


To make quantitative contact with the experiment at LENS it is crucial to take into 
account the external harmonic trapping potential V{x) = mu:‘^x^/2, which is necessary to 
prevent the fermions from e scaping . We dehne a dimensionless parameter 6 = Na\j^/a 


to describe the interactions 


ho 


l4. Il8 |. Here aho = -^/h/(men) is the characteristic oscillator 


length in the axial direction. It is somehow counterintuitive that 5 S> 1 corresponds to the 
weakly coupling limit, while 5 -C 1 corresponds to the strongly interacting regime. For a 
large number of fermions, which is about V ~ 50 experimentally j^, an efficient way to take 
the trap into account is by using the local density approximation (LDA). Together with 
the exact homogeneous equation of state of a ID multi-component Fermi gas, this gives 
an asymptotically exact results as long as A S> 1. The LDA amounts to determining the 
chemical potential /i from the local equilibrium condition 


18 


22 


231 , 


hhom [n{x)] + -mu‘^x‘^ = fig, (20) 

under the normalization restriction N = n (x) dx, where n (x) is the total linear number 
density and is nonzero inside a radius xp. We have used the subscript ‘^g” to distinguish the 
global chemical potential fig from the local chemical potential /ihom- Rewriting /Xhom into the 
dimensionless form fi[y (a:)] and 7 (x) = 2 /[n (x) aip], we hnd that 

2 ^ V [7 (a^)] + X = fig. (21) 

We solve the above LDA equations numerically jl^ . 

Fig. m reports the numerical results of the density distributions at different number of 
component k and at three typical interaction parameters Na\jy/a\f^. The linear density and 
the coordinate are in units of the peak density nTF,K = \/2Nk/ (vra/io) and the Thomas-Fermi 
radius xtf,h = \/2A /naho of an ideal gas, respectively. With increasing interaction param¬ 
eter as shown in each panel, the density distribution changes from an ideal gas distribution 
to a strongly interacting Tonks-Girardeau prohle. At the same interaction parameter, the 
density distribution become flatter and broader as the number of components k increases. 

To show explicitly the effect of high-spin bosonization, in Fig. [5] we plot again the density 
distributions of ^-component Fermi gas with varying k at a given dimensionless interaction 
parameter Na\jj/a\g = I, and compare them with the result of an interacting spinless Bose 


lO 
















FIG. 5: (Color online) Evolution of density distributions with increasing the number of components 
at the interaction parameter Na\j^/a\^ = 1. 


gas, which is obtained by replacing /Xhom with /ihom = Eq. (USD for 

24| . The profiles are shown in units of the peak density htf = \/2A^/ (Tra^o) and the 
Thonias-Fernii radius xtf = V‘2Naho, for the purpose of comparison. With increasing k, 
the profiles converge quickly to the density distribution of an interacting spinless Bose gas 
(thin line), as a result of high-spin bosonization. 


V. LOW-LYING COLLECTIVE MODES 


Experimentally, a useful way to characterize an interacting system is to measure its low- 
lying collective excitations of density oscillations j^|. Quantitative calculations of the low- 
lying collective excitations in traps can be based on the superfluid hydrodynamic description 
of the dynamics of the ID Fermi gas 2g. In such a description, the density n {x, t) and the 


velocity field v {x, t) satisfy the equation of continuity 

+^[n(3;,t)n(a;,t)] = 0, 


( 22 ) 
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and the Euler equation 


dv d 


hhom + Vtrap {x) + -TTW^ 


= 0 . 


(23) 


We consider the fluctuations of the density and the velocity held about the equilibrium 
ground state , 6n{x,t) = n{x,t) — n{x) and 6v{x,t) = v{x,t) — v{x) = v{x,t), where 
n{x) and v{x) = 0 are the equilibrium density profile and velocity held. Linearizing the 
hydrodynamic equations, one finds that |^|. 


02 


1 d 


—Sn {x,t) = -^<n— 


dt^ 


mdx \ dx 


d 


^hhom (^) 

On 


Sn {x, t) 


(24) 


The boundary condition requires that the current J{x,t) = n{x)6v {x,t) should vanish 
identically at the Thomas-Fermi radius x = ±xtf- Considering the j-th eigenmode with 
Sn {x, t) = Sn (x) exp [iujt] and removing the time-dependence, we end up with an eigenvalue 
problem, he.. 


1 d 

- r- < n 

m dx 


A 

dx 


df^hom (‘^) 
dn 


Sn ix) 


+ uj'^.Sn (x) = 0. 


( 25 ) 


We note that the above hydrodynamic description is applicable in a collisional regime 
characterized by the condition N(1 — 7~) 3> 1, where T is the transmission coefficient 


for a ID collision of two fermionic atoms along the a:-direction 16|. The calculation of the 
transmission coefficient T has been given by Olshanni in his seminal work Iq. By estimating 
a collision wavevector and by using the experimental parameters N ~ 100 and 

(^p/ci 3 D — 4.5 at LENS jsl , we find that T ~ 0.3 and A^(l — T) ~ 70 S> 1. Hence, the 
collisional regime is well reached in the recent LENS experiment |5|. 

To solve the hydrodynamic equation (|25l) . we use the powerful multi-series-expansion 


method developed in Ref. 


14| . The resulting low-lying collective mode can be classified by 


the number of nodes in its eigenfunction, he., the number index “j”. The lowest two modes 
with j = 1,2 are the dipole and breathing (compressional) modes, respectively, which can be 
excited separately by shifting the trap center or modulating the harmonic trapping frequency. 
The dipole mode is not affected by interactions according to Kohn’s theorem, and has an 
invariant frequency precisely at ui = u. Therefore, the mode frequency of the breathing 
mode provides the first means to probe the non-trivial thermodynamics of our interacting 
system. 

In Fig. ini we show the breathing mode frequency as a function of the dimensionless 
interaction parameter S = Nafj^/al^ at different number of components k. In the weak 
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FIG. 6: (Color online) Breathing mode frequency as a function of the dimensionless interaction 
parameter at different number of components. In the infinitely large number of components, the 
mode frequency will approach to that of a ID interacting spinless Bose gas (thin line). 


coupling limit {6 !:§> 1) the cloud behaves likes an ideal Fermi gas, whose breathing mode 
frequency is 2u. In the strongly interacting Tonks-Girardeau limit {S -C 1), the cloud is 
fermionized and the mode frequency is again given by 2u. Therefore, for a given number 
of components, the breathing mode frequency exhibits an interesting dip when the system 
crosses from the weak over to the strong coupling regime. We note that for k = 2, such a 
dip structure was predicted earlier by using a sum-rule approach 18|. With increasing the 
number of components k, the mode frequency decreases and hnally approaches to that of 
a ID interacting spinless Bose gas. This high-spin bosonization behavior is highlighted in 
Fig. [71 It should be noted that the Bose gas limit is very difficult to reach in the weakly 
interacting regime when <5 1. 

In Fig. [H we compare our predictions for breathing mode frequency with the experi¬ 


mental data reported by the LENS team j^. The experiment was performed at an average 
interaction parameter Na\j^/a\^ = 0.44 ± 0.08. The solid circles with error bars are the 
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FIG. 7: (Color online) Illustration of high-spin bosonization in breathing mode frequency. Here 
we show the evolution of breathing mode frequency as a function of the number of components at 
the interaction strength = 0.1 (black empty squares), 1 (red crosses), and 10 (blue solid 

circles). With increasing the number of components, the breathing mode frequency approaches to 
the frequency of a ID interacting spinless Bose gas (thin lines). 

experimental results. The empty squares with error bars are the theoretical results. The 
error bar in the theoretical result is due to the uncertainty in the interaction parameter 
^{Nalo/al^) = 0.08. The agreement between theory and experiment, within a relative 
discrepancy of a few percents, is impressive, as there is no any adjustable parameter. When 
the number of component increases, both theoretical and experimental data approach to 
the result of a ID interacting spinless Bose gas, as indicated by a thin horizontal line at 
the right part of the figure. This could be viewed as an experimental proof of the high-spin 
bosonization phenomenon. 

Qualitatively, collective modes with larger number of nodes {j > 2) exhibit the same 
feature as the breathing mode. An example is shown in Fig. [HI for the frequency of the 
third collective mode. Experimentally, however, these modes are more difficult to excite and 
measure. 
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FIG. 8: (Color online) The 3rd mode frequency as a function of the dimensionless interaction 
parameter at different number of components. In the infinitely large number of components, the 
mode frequency will approach to that of a ID interacting spinless Bose gas (thin line). 

VI. CONCLUSION 

In summary, we have investigated the thermodynamics and collective modes of ID re¬ 
pulsively interacting Fermi gases with high-spin symmetry, based on the exact Bethe ansatz 
technique beyond the mean-field method, in a homogeneous environment. This has been ex¬ 
tended to include a harmonic trap, by using the local density approximation. The equation 
of state of the system has been discussed in detail, as well as some dynamical quantities, 
including the sound velocity and low-lying collective modes. 

We have compared our collective mode prediction with a recent measurement performed 
at LENS in a Fermi gas of ^^^Yb atoms, confined in one dimension by using two-dimensional 
optical lattice. We have found excellent quantitative agreement. In addition, we have 
predicted that as the number of spin components k, increases, the mode frequency of the ID 
repulsively interacting Femi gas approaches to that of a ID interacting spinless Bose gas. 
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This intriguing high-spin bosonization phenomenon is qualitatively verified in the experiment 
in the regime with an intermediate interaction strength. 
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